Universal Critical Dynamic Form of the Vortex-Lattice Melting Line 
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A modified form of the vortex-lattice melting line is arrived at by incorporating the effects of 
critical behavior at the melting transition. Beginning with the universal form established by Blatter 
and Ivlev [Physical Review Letters 70, 2621 (1993)] which includes both thermal and quantum 
fluctuations, we then use the vortex relaxation time, r^, of a vortex-glass with a finite transition 
temperature that follows from the scaling theory of Fisher, Fisher, and Huse [Physical Review B, 
43, 130 (1991)]. This new form of the melting line is demonstrated to fit over the entire melting 
line of Yi_a;Pra;Ba2Cu306.97 (x = - 0.4) samples within the temperature range 0.03 Tc < T < Tc 
{H < 45 tesla) implying no crossover in dynamics from 3D to 2D. Generically it can be seen that a 
change in dimensionality of the vortex fluctuations along the melting line must be accompanied by 
a corresponding change in the critical exponents u and z. Such a change is observed for the highly 
anisotropic cuprate superconductor Bi2Sr2CaCu2 08. 



Since the discovery of the high-temperature layered 
cuprate superconductors, the physics of the mixed state 
of strongly type-II superconductors has received consid- 
erable attention. A significant amount of effort has been 
devoted to developing a consistent theory of the melt- 
ing transition of the vortex lattice. As there are many 
processes to be taken into consideration; thermal and 
quantum fluctuations, pinning mechanisms, correlated 
and un-correlated disorder, anisotropy, and coupling of 
the vortex lattice to the underlying electronic structure; 
arriving at an expression that is relevant over the entire 
vortex lattice melting line has proven elusive. In this 
Letter, we combine two of the most well known theo- 
ries in the physics of vortex dynamics: (1) The univer- 
sal form of the vortex lattice melting line of Blatter and 
Ivlev 0, 0i in which quantum fluctuations are taken 
into account along with the classical thermal result from 
the Lindemann criterion; and (2) the vortex glass critical 
dynamics and scaling theory of Fisher, Fisher, and Huse 
(FFH) 0, via the intrinsic vortex relaxation time, t^, 
expressed through the frequency dependent form of the 
melting line of a vortex glass with a finite melting tem- 
perature, Tg 0, 0| . The resulting expression is expected 
to be valid over the entire Hg(r) melting line, with any 
changes in the critical behavior of the vortex lattice along 
the melting transition refiected in the value of the critical 
exponent v{z + 2 — d). We use the terms vortex lattice 
and vortex glass interchangeably due to the different ter- 
minology used in 1, 2] and 3]; however, the expression 
of the melting line derived below is that of a second order 
vortex glass transition. 

We show that this new expression does in fact 
describe the entire vortex glass melting line of 
Yi_2,Pr2,Ba2Cu306.97 (x = - 0.4) thin film samples, 
YBa2Cu306.5 single crystal samples, 

(0.03 < T/Tc < 1) 0, and the electron-doped cuprate 
Smi.85Ceo.i5Cu04_j, with a constant value of the where 9 
critical exponent. However, the melting line of the 9—{'Kc\/y 



highly anisotropic superconductor Bi2Sr2CaCu2 08, 
along which the vortex lattice undergoes a well known 
3D to 2D transition at H s:! 1 kOe , exhibits a change 
in the critical exponent from that of a 3D-XY like to a 
2D-XY like value at this field. 

The vortex glass melting line, Hg(T), of the 
Yi_2;Pr2,Ba2Cu3 06.97 thin film samples (x = - 0.4) and 
ultra high purity oxygen deficient YBa2Cu306.5 single 
crystal, were investigated in magnetic fields up to 45 
tesla. The Hg(r) line was determined from electrical 
resistivity, p{H,T), measurements with H || c. Low field 
measurements, H < 9 T, were made at UCSD, whereas 
high field measurements were made at the National High 
Magnetic Field Laboratory (Tallahassee, Fl). Further ex- 
perimental details will be included in a separate publica- 
tion 7] . The Hg (T) data are shown in figure in linear 
and semi-log plots, to emphasize the quality of the fits to 
the data by the melting line equation derived below. 

In their seminal work. Blatter and Ivlev ill, 12] included 
the contribution of quantum fluctuations to the statisti- 
cal mechanics of the vortex system of a type II supercon- 
ductor. The scope of the classical formalism based upon 
the continuum elastic theory for the vortex lattice was ex- 
tended to a dynamic formalism based on the Euclidean 
action. Combining this new theoretical framework with 
the Lindemann criterion they establish a universal form 
of the melting line, ff™(T) (= Hg(T)). 

It should be noted that two slightly different expres- 
sions are arrived at in Q and 0. In their initial work 
Q , when calculating the mean squared displacement am- 
plitude, (u^), a term involving comprcssional modes is 
dropped. In the latter work y] it is not. For the first 
case, they obtain 
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reduced temperature given by 
t), Q = [Qu/in^VGl)]nTr is a 



2 



parameter measuring the relative strength of quantum 
to thermal fluctuations, Qu = dimensionless 
quantum resistance, is the Lindemann number, 
Gi = [Tc/772(0)e^3(0)]V2 is the Ginzburg number, 
is a cutoff frequency, is the scattering relaxation 
time of the quasiparticles in the vortex core given by 
the Drude formula aN = e'^riTr/m {on is the normal 
state conductivity, n is the free-carrier density, and m 
the electron mass), and d is the distance between the 
superconducting planes, is given by = min[r2p, $7^], 
where fip ~ ^le/ fJ-eTr is the kinetic cutoff frequency 
associated with the electromagnetic contribution to the 
vortex mass, and fi^ « |A is the intrinsic cutoff given 
by the gap energy due to the creation of quasiparticles 
by vortex motion. For the latter case, they find 



H„i — 
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(f + ^i + 4se/ty 



(2) 



where now 



1), S 



and 



q = 



VG- 



riTr- Either expression can be ap- 



proximated by the power-law form ^ (1 — t)" over 
temperatures ranging from Tc down to T « 0.6 Tc- By 
estimating values for Qu and ^/Gl and leaving cl and 
flTr as fitting parameters in they find a ~ 1.45, in 
close agreement with experimental values. As pointed 
out by Blatter and Ivlev, the value of the approximate 
exponent a depends on the quantum parameter Q and 
the reduced temperature 9. This then explains the ex- 
perimentally observed increase of a as the temperature 
drops below T w 0.6 0. 

Instead of putting in constant values or approxima- 
tions for the various factors of the quantum parameter 
Q, we leave them in their exact temperature and field 
dependent forms. Additionally, we allow for two intrinsic 
relaxation times; one given by the quantized level split- 
tings of the quasiparticle states within the vortex core 
Islllfll, T™''^, and the second given by the motion of the 
vortex r^, with r™'''^ as the melting transition is 

approached \vi\ . 

With (Tat cx t^°^^, then we identify r" as the relevant 
relaxation time to be associated with the cutoff frequency 
~ \/ilf./ufT^ and the expression for the quantum pa- 

ilr^. The divergent behavior of t" 



rameter Q =^ ■ 



then leads to flp < fli. (For YBa2Cu307_5 Blatter and 
Ivlev find ftp ~ 10 f2i using the assumption — t^°''^ 
The vortex relaxation time in the critical re- 
gion of the melting transition is expressed through the 
frequency dependent form of the melting line from the 
critical scaling theory of FFH. This choice is appropriate 
as we are seeking an expression relevant in the region of 
melting of the lattice. This temperature and field depen- 
dent form of Q = Q[H, t] (or q = q[H, T]), inserted back 
into equation Q (or Q), results in an equation that ties 
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FIG. 1: Vortex melting line data for 

Yi_^Pr^Ba2Cu3 06.97 films (IkOe < H < 450 kOe) 
and a YBa2Cu306.5 single crystal (100 Oe < H < 450 kOe) 
with fits of the modified melting line equation 11311 shown in 
linear (a) and semi-log (b) plots to emphasize the quality of 
the fits over the entire H and T range. 



TABLE I: Values of the Lindemann number cl, critical expo- 
nent s = v{z + 2 — d), and the quantum parameters Qo or go 
using equation 11311 ^ or equation ^ with 1121 1'^ for the data 
in figure Q 
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together the dynamics of thermal and quantum fluctua- 
tions with the critical behavior of the vortex-glass at the 
melting transition. It readily follows then that a devi- 
ation from this melting, i.e, a 3D - 2D transition, must 
be accompanied by a corresponding change in the critical 
exponent s = ^{z + 2 — d). 

Starting with the expression for the kinetic cutoff fre- 
quency, using the Bardeen-Stephen expression for the vis- 
cous drag coefficient |12j , 



(3) 



using the electromagnetic contribution of the vortex mass 
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and including the temperature dependence of 
^ — £^o/{l — t)^ , we then have 



Co 



(l-t)% (5) 



where is a shielding length, which is a few times that 
of [kpy^ 

Within the FFH critical dynamic theory of a thermo- 
dynamic vortex-glass transition [J , physical proper- 
ties are expressed through a diverging correlation length 
CvG \T—Tg\~'^ and correlation time tvg Cvg- From 
the scaling form of the frequency dependent conductiv- 
ity o-((jj) « CvG~''' S±{u!Cyq), the vortex relaxation rate 
for a glass model with a finite transition temperature ex- 
pressed as a frequency is found to be 



g ■ 



(6) 



where s = v{z + 2 — d) and d is the dimensionality of the 
vortex glass system. Returning to a relaxation time, we 
have, 



(7) 



Combining Q with (|SJl, with all expressions evaluated 
at the melting temperature Tg = T,„, we have 



4cAd I nfipaN f Tg 
Co V To \T, 



-s/2 



(l+s)/2 



(8) 



Next we use the expression for the dimensionless quan- 
tum of resistance given in jl5| . 
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FIG. 2: Fit of equation (O to the melting line, Hg{T), 
of a thin film of the electron-doped cuprate superconductor 
Smi.85Ceo.i5Cu04_a 5|. Inset is the same data on a semi-log 
plot. 

100 [ BijSr^CaCu^Oj crystal 

V T_^ = 86.0K 

"lO r "^^ ^ H^^(O) =53.5 kOe 

c =0.36 ±0.02 
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FIG. 3: Fit of equation 11311 indicating a 3D - 2D transition 
at H2D ~ 1 kOe. A corresponding change of the exponent 
v{z -\-2 ~ d) from a 3D-XY like value to that of a 2D glass 
is observed. The upper portion of the curve is fit to the 2D 
melting line expression of Schilling et al. |3 for comparison. 
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(9) ©. Combining ((U with HH), we have 



and the field dependent expression of the Ginzburg num- 
ber [li. 



(10) 



Combining |7J|, ||SJ|, ©, and (|10|l . we arrive at the final 
expression for the quantum parameters Q and q, 



Q{Hg,Tg)^^^^tHl-tY-i, (11) 



qiHg,Tg) = ^-y^-^£^tHi-ty-i 



(12) 



with Q.0 = ^ imEiL and Qo = if These can then 
be inserted back into the melting line equations ^ and 



4i/e2(0)^(l-0 



HgiQ) = 



+ ^1 + 4(gol^oTo)^c7feytni - t) 



2-4 



(13) 

It should be noted that the expression used for the 
field dependent Ginzburg number in pU|) is arrived at 
from a melting line that follows the power-law form 

H (0) — 

Hm ^ (1 — i)2, and so is strictly valid only for 

T > 0.6 Tg. In general, for a portion of the melting line 



that can be approximated by 

/ Ha 



G,{Hg) « (G,) = 



i^c2(0) 



(14) 



Also, it should be recognized that at finite frequencies 
dispersive effects lead to 0,^3 

, , ^core /I -, ,^core\ 

r]{uj) K 'S>oPs- : "2 : -2 (15) 

(1 - lajT^""-") + {loqt^'"-^) 
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with ps = 2e I'i/'l the superfluid density. Note that for 
large frequencies [toTr ^ 1) and with ps ~ [1 — t), this 
leads to the same temperature dependence of f2 found 
in © . Alternately 0, [13 , starting from ((T3|l with the 



condition wr^ ^ 1, <j-'o''> 
tJo = cjc = eB/m, with _B 



^ 1, and using ps — ~en, 

^ Hc2 



$2 



2Tre 



(16) 



which is identical to the Bardeen-Stephen result Q. 

The melting line data, Hg(T), obtained in this study 
of Yi_a:Pra:Ba2Cu3 06.97 films and a YBa2Cu3 06.5 sin- 
gle crystal were fitted by equation l(T^ (figure (Q- a, b)) 
and that obtained by substitution of H12|l into (not 
shown). With the experimental value Qq = 0.34 ± 0.15 
obtained from equation (|13|l . we solve for the value of 
OoT^'o for the YBazCuaOy-d" film. Using = I2A, 
e = 1/8, PAT « 2 X lO'^Qcm, A « 14 meV H^, and 
with Xd a few times (1 - 3 ) kp^; ^oTro ~ 0.65. If instead 
we use the second modified melting line equation ^ with 
Pth = 2.5, then we get an equally good fit indistinguish- 
able from (|13|) with similar values of i'{z — 1) and c|, 
with Qq w 0.1. In general, it is seen (see tabled that the 
value of is reduced from Qq by a factor of ~ 1/2 - 1/10. 
The values from both fits are listed in tabled The result, 
q < Q, can be understood as a relative reduction of the 
quantum contribution to displacements by the inclusion 
of the thermal contribution of the compressional modes 
to the mean squared displacement amplitude of the vor- 
tices, (u^), used in the melting line equation which 
are not in 

Additionally, equation (|13|l was used to analyze melt- 
ing line data from a Smi.g5Ceo.i5Cu04_j, film la 
(figureHJ, and a Bi2Sr2CaCu208 single crystal |g| 
(figure 01 . The data for all of the above samples are 
described well by equation ((T!^ over the entire range 
of field/ temperatures examined, with the exception of 
the Bi2Sr2CaCu208 data, which can be fit by equation 
(|13|l in two segments, with a corresponding change in 
the critical exponent s. This latter case is interpreted 
readily as evidence for a 3D-2D vortex glass transition at 
H2D « 1 kOe, which is well established in this highly 
anisotropic compound. Since the data in the other cases 
considered here can be fit over the entire (available) range 
of the Hg (T) line with one set of critical exponents corre- 
sponding to 3D, this implies that the vortex glass remains 
in a 3D state along the entire line. 

While it is a priori unclear as to which equation for the 
melting line is more appropriate for a given system, accu- 
rate measurements of the many physical properties (Hc2, 
Gi, ^0, etc.) that enter these equations would, in prin- 
ciple, make it possible to distinguish between the two. 
Either way, we have arrived at a form of the melting line 
that unifies the quantum/thermal nature of vortex fluc- 



tuations with the critical dynamic behavior of vortices 
in the region of the melting line transition, providing a 
more complete picture of the physics involved. 
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